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1. Introduction 

From long time ago, many papers have appeared on the nonperturbative computational 
technique of quantum systems based on the idea of interpolating actions yy. The 
technique is called under various names such as variational purterbation, Gaussian 
approximation, linear delta expansion, interpolation method and so on. All of these 
methods exploit artificially incorporated parameter 6 which interpolates non-trivial 
action of interest (at 5=1) and some solvable one (at 5 = 0). When 5 is set to 
unity, the total action becomes independent of the solvable one and the independence is 
often utilized to obtain non-trivial results by the use of principle of minimum sensitivity 

In some cases of those approaches, the introduced parameter 5 can be viewed as 
the dilation parameter for one of parameters originally included. On the point of view 
of dilation, one can use the method in new ways. For example, in lattice field theories, 
we can start with the strong coupling expansion which is often valid at a large lattice 
spacing a and then perform the expansion in the dilation parameter 5 after the shift of 
a a{l — 5). This technique, simply called delta expansion, was proposed and applied 
to some lattice models and it was shown that the continuum scaling emerges in the 
strong coupling series [3l H] . The non-perturbative mass gap in the continuum limit was 
then evaluated to good accuracy. To be self-contained, we first describe the method 
below. 

Let a quantity of interest VL has an expansion: 

^ = E (1) 

Here, M denotes the mass square in momentum space rescaled to be dimensionless. 
The continuum scaling of Q is given by the behavior as M — > 0. To access the 
small M behavior of Q, we dilate Q{M) around the continuum limit by shifting M 
as M — > M(l — 6) where < 5 < 1. Suppose then 

n{M) ~ AM"", {a > 0), (2) 

in the scaling region. Dilation leads that f2(M(l — 6)) ~ AM~"(1 — and as 5 ^ 1 
dilated Q diverges for any fixed M > as in the manner that Q{M) diverges as M — 0. 
Since the original scaling behavior is thus transferred to the behavior of dilated function 
in the 6^1 limit, there is a possibility that the scaling behavior could emerge in dilated 
large M expansion of Q, even if we have only finite series to order . To put this 
idea into practice, a critical step is to expand Q{M{1 — 6)) in dilation parameter 6 both 
at large and small M in an appropriate manner. Now, i^'th order Q{M{1 — 6)) at large 
M reads f2(M(l — 5)) = Y.k=i m''(i-5)^ ■ [31 H]' the expansion in 5 is truncated along 
with the conventional rule of partial sum, 

fi((l-5)M)~ ^ijM-'5f (3) 

l<i+j<K 

Here uoij = uji,u!2j = {] + l)co'2,t^3j = -^^^^t^^^^cjs, • ■ -. The order of term M~^5^ is 
assigned as i + j and it is included as long as i + j < K. Note that maximum order of 
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5 is K — 1. On the other hand, at small M, we need only the leading term ([2]) giving 
asymptotic scaling. We expand it in 5 to 5^~^ to keep accordance with the highest 
power of (5 of ([3]). The result then reads 

m - 6)M) ~ AM- E ^^^±I):_^^±1^5^ (4) 

i=0 ^• 

Setting 6 = 1 which means the dilation to infinite amount, the truncated series of Q at 
M ^ 1 and M -C 1 were compared with each other and was found to exhibit the 
correct "scahng behavior" (jl]) at 5 = 1. We emphasize that, due to dilation, the scaling 
is observable in a wide region of shifted M of order M ~ 0(1) or even larger. 

As explicit in ([3]) the series at large M has two expansion parameters, and 
S. If both parameters are small enough, the truncation of ([3]) would be the most 
reasonable. However, 6 is always set to 1 eventually. Then, terms of different j for 
fixed i, bij{M~^y6f mix among themselves. Hence at 5 = 1 the order assignment loses 
the basis. This motivates us to seek for and consider other truncation schemes admitting 
analytic control. Within those schemes, we like to examine convergence properties as the 
order of expansion increases. In addition we will point out that when M and the order 
N of the expansion in S is large enough with M/N kept fixed, large M series approaches 
to the Borel transform of the original series Thus, by exploring various possible 

schemes, we can obtain computational flexibility toward models complex and not fully 
understood yet. As a demonstration of important role of Borel transform limit, we 
reinvestigate the continuum limit of the Gross-Neveu model with Wilson fermion. 

This paper is organized as follows: In the next section, we introduce a class of new 
truncation schemes which we call " square" schemes. Then we formally show that those 
schemes lead to Borel transform in a certain limit. In section 3, delta expansion in 
square schemes is examined in detail by applying the schemes to a simple model, the M 
component anharmonic oscillator in the large J\f limit. Then, we apply Borel transform 
to study the scaling properties and evaluate the mass gap. In section 4, continuum 
scaling of three models, ordinary anharmonic oscillator which corresponds to A/" = 1 
case, two dimensional (2d) non-linear a and Gross-Neveu model with Wilson fermion 
are studied at large A/" in the Borel transform approach. In the non-linear a model, 
to improve the accuracy of approximating the continuum limit, Symanzik's improved 
action [8j will be discussed to the first order. In Gross-Neveu model, we point out that 
the conventional truncation scheme fails to capture the continuum scaling but the new 
one works good. The results of this work is summarized in the last section. In the study 
of anharmonic oscillators, we will confine ourselves with the pure anharmonic case where 
the harmonic mass term is absent from the action. 

2. Delta expansion in square truncation and its certain limit leading to 
Borel transform 

In this section we define square schemes and study the limit where the order of expansion 
in 6 is taken to large enough. 
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To define a square scheme, let f2(M(l — 5)) be first expanded in 6 to 5^ such that 
VL{M{1 - 5)) = fio(M) + ni(M)5 + ^2(^)52 + . . . + n^(M)5^ (^oiM) = n{M)). At 
M <^ 1, expanding Qk in ^ and collecting leading terms for each k, the result (jll) with 
the replacement — 1 ^ is derived. Then, we find in the S ^ 1 limit 

Q{M{l-6)) ^ AM-'^ZNia) (5) 

where 

r(Ar + a) 
= (Ar-l)!r(« + l)- 

At M > 1, we expand fi^ in powers of to (M"^)^. Then at M > 1, the partial 
sum of dilated Q is formally written as 

fi((l-5)M)~ J2 u^ijM-'Sf (7) 

l<i<K,l<j<N 

The order of the truncated double series ([7]) are labeled by K and A^. Finally we set 
5 = 1 in ([7]). To obtain ([7j) at 5 = 1 from the original 1/M series ([1]), we suffice to expand 

dilated term M-'=(l - 5)-*^ to the order A^ such that M-'=(l + + k{k + 1) /2!(52 H h 

+ 1) ■ ■ ■ (A; + N)/N\6^). Then, setting 5 = 1, we obtain the simple transformation 
rule, 

M-'^M-'{^^^). (8) 
This gives the dilated f2 at 5 = 1 to orders {K, N) as 

^iM'^ k ^ ' ^^iM'^J' 

We always deal with the truncated series and there appears no divergence connected 
to the M limit even when 6 is set to unity. Though fi(0 x M) is M independent, the 
right hand side of ([5]) and have M dependence. As the order of expansion increases, 
the residual M dependence should become weaker. 

In 2d non-linear a model, we consider the scaling of bare coupling {Q = -). Then 
a = in ([5]) and the scaling behavior is a logarithmic one, Q{M) ~ AlogM + B. Since 
logM turns to logM(l — 5) under the dilation and is expanded as logM — J2k=i T' 
have 

N ^ 

logM ^ logM (10) 
fc=i 

where S has set to 1. Hence, we have the following modification of the continuum 
logarithmic scaling, 

N -1 

AlogM + B ^ A(^\ogM -J2t) + B. (11) 

k=i ^ 

By calculating and comparing it with ([5]) or ( ITTi) . the informations on A and B can 
be extracted. 



Continuum scaling in expansions effective at a large lattice spacing 



5 



Now, we consider the limit that the order is taken large enough. In the limit, 

+ 
k 



the lactor ( 1 behaves as 



Hence M"*^ in original expansion transforms to ^{M / N)~'' . We then obtain 



K ^ K 



If M is 0(1), all terms in the right-hand-side (RHS) of (fT3|) diverge. However, if M 
is taken large as the same order of A^, the series in the RHS of (fT3l) becomes sensible. 
Thus, it is tempting to consider the following correlated limit, 

M 

M,N oo with — fixed. (14) 

By defining a new variable M by 

M := f , (15) 
we summarize the result of square scheme to large orders as 

It is obvious that the obtained result is nothing but that of Borel transform with respect 
to M. On the other hand, at small M, we find 

niOxM)^Ax M- \,^\ . (17) 

i (1 + a) 

in the limit ( |T4|) . For the case of the logarithmic scaling, we find from Y^k=i \ ~^ 
\ogN + 7£; {pfE stands for the Euler constant) that 

logM logM-7i5. (18) 

Hence, 

Vlif) X M) A{\og M -^e) + B. (19) 

The above argument leading to Borel transform is a formal one and needs detailed 
explanation. Taking the M component anharmonic oscillator as an example, we like to 
demonstrate that the limit (fT^ is legitimate and sufficient for our purpose (see the next 
section) . 

Note that is completely combined with the mass parameter M. The effective 
variable is M = M/N and we need not specify the value of iV, which plays the role 
of the regulator of the delta expansion. By the comparison of ( fT6i) to ( fT7|) or ( fT9i) . we 
can obtain the informations of the exponent a, constants A and B, if ( IT6|) exhibits the 
scaling behavior at a region of small M. 
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3. Application to the M component anharmonic oscillator 

As the first application of square delta expansion and Borel transform (BT) to concrete 
models, we address to M component anharmonic oscillator in the large limit, since 
the model is solvable within self-consistent approach. The lattice spacing is a and 
M component field (pn = (0^'^, 0^, ■ ■ ■ , 0^) denotes the dynamical degree on a cite n 
{n = 0, ±1 ± 2, ■ ■ ■ , ±Z). The massless action is then given by 

C \^ rl/0„+l-0n\2 A .7212" 



(20) 



n=—L 

The action can be rewritten by rescaling fields from (p to = (aA/4)^/^0, giving 



S = E_(M + j^i^lV) -PT._^n- (21) 



-L n=—L 



where 



The second term represents the hopping term by which nearest neighbour fields tend to 
align 3. 

Now, the self- consistent method gives 

2 > dp 1 

m = X f 5— r-, (23) 

j-vr/a 27r ^(1 — COS pa) + 

where m stands for the dynamical mass at spacing a. Changing the integration variable 
from p to 6 hj p = 6/a and defining the dimensionless mass variable by 

M := (ma)^ (24) 

we find 



dd 1 



1/2 

(25) 



M J-n 27r2(l -COS0) +M 
It is easy to obtain exact /5 as a function of M by calculating the integration explicitly 
but it is not our aim. The starting point of our argument must be of generally available. 
It is a series expansion of P in powers of 1/M which is obtainable from (1211) . However, 
in the present case, the result can be obtained easier from fl25l) : 

^ M 4M5 ti ' 

We emphasize that the self consistent method is used just to produce the large M 
expansion in an efficient way. In the next subsection, using only the above series, we 
can address to the scaling behavior, 

(3 ~ V2M-^/\ (27) 

which is derived also from (l25l) . 

I The expansion parameter /3 in our hopping expansion is twice of k which is used as the conventional 
hopping expansion. 
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3.1. Delta expansion in a square truncation 

In this subsection, we use the method of delta expansion in a square truncation scheme 
to capture the scahng behavior ( 1271) in the large M series. We proceed by correlating 
the orders of 1/M and 5 expansions by 

K = N. (28) 

Then, to the full order K., we have 

m-ty't')^.. (29) 

fe=l 

Figure 1 shows the plots of [3 fl2Bl) and D[ff\ fl2Ul) from 2nd to 8th orders. As explicit 
from Figure 1(2), D[(f\ shows rough scaling at rather large M. 




Figure 1. (1) Plots of (i from 2nd to 8th orders. (2) Plots of D[/3]/Zk(3/4) from 
2nd to 8th orders. In both figures, each dotted line represents the asymptotic scaling 
/3 ~ y2M-3/4. 

To access the continuum scaling, it is more convenient to deal with the logarithm 
of (3. This is because when the small M behavior of (3 is such as (3 = y4M~^/^(l + ciM + 
C2M^ + ■ ■ ■) 5 delta expansion leaves all corrections to the asymptotic scaling to any 
finite orders in 5. However, for log 13 = log A — 3/4 log M + ciM + O(M^), the low order 
corrections, M^, M^, ■ ■ -, are dilated and expanded as M(l - 6), M^(l - 26 + 6"^), ■ ■ ■, 
and they disappear by setting 5 = 1 at several orders. Now log P at large M becomes 
2 1 2 16 16 256 2 b' 

and 

nri m i ^ ^ ^ / + 1 x 1 / + 2 x 2 / K + 3 x 16 
^[log/3]=log- + E--( ^ )- + ( 2 )jf^-{ 3 )^ + --- 

It should be reminded that if (3 is of order K, then the corresponding log (3 becomes of 
order K — 1. 
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At small M, log /3 ~ log ^2 - | log M and D\[og (3\ ~ log ^2 - |(log M - ^f^^ i). 
The behavior of D[log/5] is plotted in Figure 2. It is clearly shown that we can observe 
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Figure 2. Plots of exp(i:)[log/3] - fEfcii i) from 1st to 7th orders. The dotted hne 
represents the asymptotic scaling. 

the scaling behavior in 1/M series more explicitly. 

To confirm the value of the exponent a, we deal with This function behaves 

at small M as 

d\og[3 

— = —« + ••• 32 

(91ogM ^ ' 

where ■ ■ ■ represents the terms which vanish in the M — > limit. At large M it behaves 
as follows: 

aiog/3 _ _ 1 4 

91^ " ~^ + M " + ■ ■ ■ • ^^^^ 
If we introduce 5 to dilate the region of M by M ^ M(l — 5) and take 5^1 
limit, -P[ g/°g jfj ] approaches to the uniform function taking the value —a. We note, 
however, when the dilated function is obtained as an approximant, it would have weak 
M dependence and a plateau suggesting —a. By performing a square delta expansion 
on the large M series of , we have obtained the behavior shown in Figure 3. At 

several orders, the truncated large M series exhibits almost uniform behavior at large 
M, and the values at the plateaus are close to the exact value of —a = —3/4. Thus, 
the square delta expansion indeed shows correct value of a. 

The amplitude a/2 in (1271) can be evaluated as well. From (1271) . log/? + |logM = 
Q{M) behaves in the scaling region, 

g(M) = log V2 + ---, (34) 

where ■ ■ ■ stands for the corrections all of which vanish in the M ^ limit. The dilation 
around M = makes the behavior of Q{M{1 — 6)) stationary within a wide region at S 
close to 1 and, as 5 — 1, Q{M(1 — 5)) converges to log -\/2 at any finite M. 

Then, we like to show that the delta expansion on the large M series of Q recovers 
the flatness and emerged plateau indicates the correct value, log ^/2. Figure 4 shows the 
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Figure 3. Plots oi D[§^^] at 2nd, 5th, 8th, 11th, 14th and 17th orders. The dotted 
hne represents —a = — |. 

plot of 

DIQ^ ^ .og2 - l(logM - gj) - ( 1 )-L . ( '^■^^ ^ )^ . . . . (35) 

at 2nd, 5th, 8th, 11th, 14th and 17th orders. At odd orders there exists only one 
extremum value and it may be regarded as the representative values on the plateau. 
Then, as proposed in [2], it is natural to take the extremum value as the approximation 
of log A = log a/2 = 0.346574 ■ ■ -. The results at 5th, 11th and 17th orders are 

log A: 0.3098, 0.3397, 0.3450 (36) 

and these values occur at M = 13.4558, 15.7481, 16.7, respectively. Our results above 
several orders are found to be in good agreement with the exact value. The mass gap 
can be computed by using above approximants for log A by m = (y4/2)^/^A^/^, though 
we omit explicit results. 



. 8 
. 6 

D[Q] 

. 4 
. 2 



5 10 15 20 25 30 

M 

Figure 4. Plots oi D[Q] = D[\og (3 + 3 / Alog M] at 2nd, 5th, 8th, 11th, 14th and 17th 
orders. The dotted line represents log — 0.346574 • • •. 
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Table 1. Approximation of the scaling amplitude \og\/2 ~ 0.346574 •• -in various 
truncations in the square delta expansion. The value of M at which Q reaches 
extrcmum is written in the parenthesis. 



order K 


N = K 


N = K+l 


N = K + 2 


N = K + 2, 


5 

11 

17 


0.30985 (13.456) 
0.33974 (15.748) 
0.34501 (16.700) 


0.31341 (15.305) 
0.34029 (16.783) 
0.34513 (17.420) 


0.31601 (17.152) 
0.34075 (17.816) 
0.34523 (18.140) 


0.31799 (18.996) 
0.34113 (18.848) 
0.34532 (18.859) 




order K 


N = K + A 


N = K + b 


N = K+IQ 


N = K + 20 


5 

11 

17 


0.31954 (20.839) 
0.34145 (19.879) 
0.34539 (19.578) 


0.32079 (22.680) 
0.34173 (20.910) 
0.34546 (20.296) 


0.32457 (31.874) 
0.34267 (26.056) 
0.34569 (23.883) 


0.32760 (50.239) 
0.34356 (36.328) 
0.34594 (31.041) 



order K N = K + m iV = is: + 40 

5 0.32891 (68.593) 0.32963 (86.943) 

11 0.34397 (46.587) 0.34421 (56.840) 

17 0.34606 (38.189) 0.34613 (45.331) 



3.2. Other truncations in square schemes 

In the previous subsection, a square truncation was investigated under the condition 
N = K. There are other infinite choices of the truncation rule in the class of square 
schemes. We here study the results of choices that 

N = K + L, (37) 

where L will be assigned some positive integer. We write here the result of evaluating 
the amplitude of scaling, log \/2, which is the leading term of Q at small M. The results 
are summarized in Table 1. 

We find that as L increases the approximant increases monotonically. Further, 
approximation is improved as L increases. Note that the value of M at which the 
function D[Q] becomes extremum also grows with N. Let the value of M at which 
D[Q] takes extremum value be M*. Then, at = 17, the ratio M* /N has the following 
values at L = 0, 1, 2, ■ ■ ■ , 5, 10, 20, 30 and 40, 

0.982355, 0.967794, 0.954736, 0.942958, 0.932281, 

0.922556, 0.884561, 0.838967, 0.812543, 0.795287. (38) 

The ratio is of order 0(1) and decreases gradually. Now, it is apparent that, in small M 
region, the delta expanded 1/M series cannot offer us the scaling. Rather, scaling can 
be observed at the region of M ~ 0{N). Hence the small M region can be neglected 
and the correlated limit turns out to be a natural limit to consider. This is the 
reason behind considering BT limit. 
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We emphasize that scahng reveals itself at M ~ 0{N) is a characteristic feature of 
square schemes. In the conventional scheme [3lll], the scaling is observed at M ~ 0(1). 

3.3. Results in Borel transform limit 

In this subsection, we study the scaling behavior of (3 via 1/M expansion by the use of 
Borel transform. 

The Borel transform of log (3 reads at large M 



log ^+1E 

M 



1 

M 



2!M2 



2 °° 6' 

+ ■ ■ ■ = log ^ + 7b + V — ^ 



(39) 

and 8th. It is 



The behavior of fl5^ is plotted in Figure 5 at orders 2nd, 3rd, ■ 
obvious from (139|) that Borel transform improves the small M behavior: Comparing 
(!39l) with( |30l) . i?[log/3] has the inverse factorial in the coefficient and its convergence 
radius is larger. In the model under consideration, log (3 converges for < < i. 
Hence i?[log/5] has infinite radius of convergence. Thus the scaling behavior may be 
seen in 1/M series. 



B[log p] 



2ndMV 

8th\Vv 




3rd/// 





. 5 

A 

log M 



Figure 5. Plots of B[log/3] from 2nd to 8th orders. 
log[\/2Af-3/Vr(3/4-t-l)]. 



The dotted hne represents 



The evaluation of the exponent a is easily done by considering the Borel transform 
of qI^^^j . We here omit the task and directly turn to the evaluation of the amplitude 
A which is connected to the magnitude of mass gap. We consider Q = log (3 + | logM 
and its Borel transform. Figure 6 shows the plots of B[Q] = log 2 — ^(logM — •je) + 
En=i and B[Q]\^^, = log V2. 

In terms of M, the transformed series does not possess a plateau of wide region. 
This is a natural result since values of M of 0{N) correspond to values of M of 0(1). 
In other words a wide plateau at M ~ 0{N) in square schemes are mapped to a 
region of M ~ 0(1). Since the stationarity is conserved by the mapping, however, we 
approximate log A by stationary value of B[Q] . The results at orders 5th, 11th and 
17th follows: 



log A: 0.332143, 0.345087, 0.346381. 



(40) 




Figure 6. Plots of B[Q] from 2nd to 8th orders. The dotted hne represents logv^- 

These values are realized at M = 1.834, 1.024, 0.712, respectively. It is now sure that the 
above sequence approaches to logv^ = 0.346574- ■ -, the exact value of the amplitude. 
The corresponding mass gap is computed as m/X^^^ = 0.7861, 0.79291, 0.79359 and in 
good agreement with the exact value, m/\^^^ = 2^^^^ = 0.793701- ■ ■. Comparison of 
the above approximants to those in square schemes at finite L clarifies that BT limit 
produces best results. 

4. Application to other models 

Also in models we are going to discuss from now on, BT limit in a class of square schemes 
gives most accurate result than any of finite N. The reason would be that the power 
like correction {k > 0) remains to A; — 1th order square scheme but, in BT limit, 
it vanishes at every order. This means that the correction to the asymptotic scaling is 
smaller in BT limit. Thus, we report the results only in BT limit. 

4.1. Anharmonic oscillator 

In this subsection we apply the delta expansion in the BT limit to the single component 
anharmonic oscillator. At A/" = 1, the anharmonic oscillator cannot be solved and only 
numerical results are known to high accuracy. It therefore serves us a good theoretical 
laboratory to examine our method. 

In terms of the rescaled field on a site n, the action S reads 

5 = E ^iVn) - /5 E ^n+l¥>n, ^(^) = V + ■ (41) 

n n 

Even when the system is massless at the level of action, the fields (po and ipn at large 
separation has the finite correlation length. Then, Fourier transform of the two point 
function {(foipn) at large n defines M corresponding to IHM . By calculating {ipoipn) by 
the use of hopping expansion [9], one can obtain M as a series in (3. Then, inverting 
and (3, we have large mass expansion of /3. In [4], is obtained up to M~^. 
Unfortunately, as in the large N case, the Borel transform of (3{M) does not show clear 
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sign of the continuum scaling in its effective region even at 9th order. However, taking 
the logarithm of 13 improves the state of affairs. We find 

log/3 = -log(pM) + (-3 + -^)-^ + (6--^+ ^ ' ^ 



- . 



4p2yM V 4p2 I2p4y]vf2 

33 47 27 5 x 1 /411 431 983 457 

+ 



/ 33 47 27 5x1 / 



2 8p2 32p4 128p6/M3 V 8 16p2 leOp^ 768p6 
385 X 1 / 3417 973 1553 703 1303 61 



18432p8/M4 V 20 8p2 40p4 120p6 3072p8 5120p 
1 /2367 2173 72411 45199 47081 1391 

XTT7+ ^ + 
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M5 V 4 4p2 320p4 9g0p6 8960p8 4608pio 

583 X 1 / 117981 38511 800529 85869 3667987 



^82944pi2/M6 ^ V 56 ^ 16p2 640p4 ^ 256p6 71680p8 

387119 62843 52195 x 1 /489681 678201 

~86016pio ^ 294912pi2 ^ 12386304^^ ^ ^~64 64p2 

3409315 11285453 1387540351 7581293 10207705 

+^T7r-. .... . + 



512p4 5120p6 3225600p8 147456pio 2752512pi2 
350285 861575 x 1 q. 



2359296pi4 339738624pi6/ 

h' 



oo t/ 

-log(pM) + E77^, (42) 

k=l 



where p = r(3/4)/r(l/4) = 0.337989 • • •. Then at large M, Borel transform of log/3 is 
given by 

°° h' ^ 

B[\ogP] = - logp - logM + + E ttM-^ (43) 

k=i ■ 

Here, truncation of -B[log/?] up to will be called Kth. order approximant. Now, 

the scaling behavior of log 13 is given by 

3 

log ^ ~ log A - - log M, (44) 
and its Borel transform reads 

B[\ogP] ~ log A - ^(log M - 7£;), (45) 



where log A is known numerically as log A = 0.81841 ■ ■ 

Figure 7 shows the plots of i?[log/3] at 1st and 8th orders. At 8th order, we find 
that the scaling behavior is observed around logM ~ in the Borel transformed 1/M 
series. Having captured the scaling behavior, we can evaluate the constant A which 
leads to the mass gap in the continuum limit by m = (A/2)2/^A^/'^. The evaluation 
step goes as in the previous section. We deal with Q = log /? + | log M and consider its 
Borel transform. B[Q] is plotted in Figure 8. We find that the effective region grows 
to the smaller M region as the order increases. The value of log A is indicated by the 

§ The value of log A is calculated from the value of the mass gap referred in [lOj 
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Figure 7. i3[log/3] at 1st and 8th orders. The dotted hne represents log A— |(logM — 
7_e) where log A = 0.81841 • • •. 

stationary value as in the case of large M limit. At i^T = 1, 3, 5, 7, we have following 
approximants of log A: 

log A: 0.68467, 0.75734, 0.78197, 0.79358. (46) 

Since the exact value of log A is 0.81841 ■■ ■ we can say that delta expansion in BT 
limit is successfully working. The mass gap m is then computed at respective orders as 
follows: 

TfX 

—r : 0.9943, 1.0437, 1.0610, 1.0692. (47) 

The exact value of m is known to be m = 1.087096 ■ ■ ■ x X^^^ [10] and the results at 
several orders give good approximation. 

1 . 4 
1 . 2 

B[Q] 

. 8 
. 6 
. 4 
. 2 

2 4 ^6 8 10 

M 

Figure 8. Plots of B[Q] = B[logf3 + flogAf] at K ^ 1, ■ ■ ■ ,8. The dotted line 
represents the value, log A = 0.81841 •• •. 

4-2. 2d non-linear a model 

Up to now, we considered quantum mechanical cases where the models are defined on 
one-dimensional lattice. In this subsection, we discuss a model field theory, the non- 
linear a model in the large Af limit at two dimension. The 2d non-linear a model can be 
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solved out in the large Af limit and enjoys interesting properties of asymptotic freedom 
and dynamical mass generation. We study the continuum scaling behavior via the large 
mass expansion under BT limit. 

The 2d non-linear a model on continuum Euclidean space is defined by the 
action, C = jj {d^aY where / denotes the bare coupling constant and the fields 
a^{x) {A = 1,2,---,M) obey the constraint, a^{x) = Ea=i = A/". The 

discretized space we work with is the periodic square lattice with the lattice spacing a 
where a site is labeled by two integers, (ni,n2) = n. A simple version of the lattice 
action is given by 



where 



S = 2PY.^l-Pll E ^n-C?„+e„ (48) 
n n fj.=l,2 

P--=j (49) 

and we call P as the hopping parameter. In ( l48l) . (Tn+e^ stands for the nearest neighbour 
spin of (Jn with ei = (1,0) and e2 = (0, 1). The constraint is written as c?^ = A/" and 
the first term in (B8|) is actually a constant that can be omitted. 

In the large A/" limit, the model can be solved out and f3 is given from the constraint 

as 

(Pp 1 

^ = I. (27r)2 M + 2 E^=i,2(l - cosp^) ' 

where 

M = {may. (51) 

We start our analysis with 1/M expansion of (3, which is equally available in other lattice 
models. It is easy from (!50|) to obtain 

1 4 20 112 676 4304 7, 

Previous two models have a common feature that the logarithm of (3 is more 
convenient to address the continuum scaling than f3 itself. It does not apply to 
the present model since /3 itself behaves logarithmically in the scaling region, f3 ~ 
-^log(M/32). log/3 behaves like ~ log(-^ log[M/32]) at small enough M, and to 
capture such a behavior is more difficult. Hence, we consider Borel transform of jS 
which is easily obtained as B[i3] = Y^^Li aX M ^ 1 and at the scaling region, 

B[(5] ^(log M-^E- log 32) (M < 1). Unfortunately, the scaling behavior of B[(3] 

is seen at M ^ 1 only roughly, though remarkable improvement is found as shown in 
Figure 9(2). To improve the status, we employ Symanzik improvement of the lattice 
action [8]. The reason why Symanzik's modification of the lattice action helps us to 
capture the asymptotic scaling is described in detail in [3]. To say it briefiy, the leading 
correction to the asymptotic scaling is MlogM and this affects the small M behavior 
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Figure 9. (1) Plots of /3 in the large M expansion at 2nd, 5tli, 8tli, lltli, 14th, 17th 
and 20th orders. (2) Plots of B[(3] in the large M expansion at 2nd, 5th, 8th, 11th, 
14th, 17th and 20th orders. (3) Plots of B[(3] in the large M expansion at 2nd, 5th, 
8th, 11th, 14th, 17th and 20th orders at the first order Symanzik program. 



even at small M. Then, the term can be subtracted by introducing the next-to-the- 
nearest neighbour coupling term J2fj. ■ o^n+2e^ into the action. At the first order of 
Symanzik imporovement program, the action becomes 

S = PY. -Tn+e^ + ^ II " '?n+2e^ • (53) 

The constraint relation at large A/" now reads 

^ " (27r)2 M + 5 - I Em cosp^ + I E;. cos 2p^ ' 
In the improved action, the large M expansion and the value of the constant part in the 
scaling behavior becomes as follows: Large M expansion reads [3], 

1 5 1157 8419 ^ bk 



^ M M^^ 36M3 36M4 M'^ ^^^^ 

and the scaling behavior. 



/?~-^logM + S', (56) 

where B' = B - 0.0471699 ■ ■ ■ = 0.2286245 ■ ■ -. By plotting B[(3] in 1/M expansion, 
we find clear logarithmic scaling in accordance with the asymptotic freedom (see Figure 
9(3)). 
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Table 2. Approximants of B' in Borel transform (BT) limit. Exact value of B' is 
B' = 0.2286245 • • •. The value of M at which B[Q] becomes extremum is shown in the 
parenthesis. 



order 


B' approximants 


B' approximants 


2 
3 
4 


0.215929 (7.836) 




5 


0.213059 (6.372) 




6 


0.212669 (5.915) 


0.213555 (3.961) 


7 


0.212765 (6.066) 




8 


0.212748 (6.038) 


0.216666 (2.966) 


9 


0.212751 (6.043) 




10 


0.212751 (6.042) 


0.219193 (2.436) 


20 


0.212751 (6.043) 


0.225305 (1.340) 



Having observed the scaling, we turn to the evaluation of the constant B' which 
is directly connected to the mass gap hy m — exp(27rS')A (A stands for the scale 
parameter in the first order Symanzik model). For the purpose it is convenient to deal 
with g = + ^ log M = ^ logM + J2T=i Mk which tends to B' in the M ^ hmit. 
By Borel transforming Q we have plotted the resulting functions at 2nd, 5th, 8th, 11th, 
14th and 17th orders (see Figure 10). Approximation of B' by extremum values are 

0.5 
0.4 
0.3 

B[Q] 

0.2 

. 1 

2 4 6 8 10 

M 

Figure 10. Plots of B[Q] = B[l3 + ^ logM] at if = 2, 5, 8, 11, 14, 17.20. The dotted 
line represents the value of B' = 0.2286245 • • •. 

shown in Table 2. In BT limit, there appear two extrema at 6th order. This is a signal 
of appearance of new family of extrema at higher orders. First family of extrema stays 
at M ~ 6 to all higher orders. The second family appears from 6th order and the 
location moves to smaller M as the order increases. It is obvious that we should trace 
the second family, since it signals lira j^^_^qB[Q]. 

To summarize, the sequence of extrema at smallest M suggests strongly that it 
converges to the exact value of B' and the dynamical mass can be calculated from the 
large M expansion of /3. 
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4-3. 2d Gross-Neveu model with Wilson fermion 

The models considered up to now allow both summation schemes, conventional and 
square ones to us for simulating continuum limit. It would be nice if the truncation 
presented in this work has an essential importance for the purpose. 

Here we show a theoretical example, 2d Gross-Neveu model in the large A/" limit [6] 
with Wilson fermion where conventional scheme fails but square and BT schemes 
succeed to capture the continuum scaling. 

Gross-Neveu model describes quartic interactions between massless A/" component 
fermion fields ipa = 2, ■ ■ ■ ,J\f; a = 1, 2). Introducing the auxihary fields a on each 
site, the action on the lattice reads [H] 



^ = - -7/.)V'(n + e^) +^(n + e^)(r + 7^)V^(n)] 

+ 5:(2ar + aV(n))V^(n)^(n) + ^ Y.i^{n) - 5m)\ (57) 

where 7 matrices are given by Pauli matrices, 

7i = 0-2, 72 = 75 = 0^3 = ^7i72- (58) 
The parameter r is called Wilson parameter and as long as r 7^ 0, the 75 symmetry 
related to the transformation —>■ j^ip is explicitly broken. Thus the mass protection 
mechanism is absent and there exists mass divergence even when the current mass is 
zero. The parameter 6m represents the contribution of the mass counter term to cancel 
out the divergence. Now, to fix the mass counter term, we use the perturbative analysis 
and compute the a tad pole in the large J\f limit. We employ the renormalization 
condition that the tad pole exactly vanishes. Then we find 

5m = -^CAr) (59) 
a 

a(r)= r '^'^ ^E,(l-cosp^) 

i-.(27r)2{rE,(l-cosp,)P + E,sin2p/ ^ ^ 

Under this prescription, the mass cannot be created to all perturbative orders. Hence, 
the dynamical mass can be generated only non-perturbatively. 

One can discuss the recovery of the broken 75 symmetry in the continuum limit by 
studying the effective potential of a. In the large M limit, it is exactly calculated to 
give [II] 

= [IIsin2p^+(aL+r^(l-cosp^))2J,(61) 

where ctl = aa, 5mL = 5m ■ a and Vl = Va^ ■ The restoration of 75 sjmametry in a 
limit is confirmed by expanding Vj, in ctl. The result is 

= -(^ + 2C,)a. + - ^0 + 2C.)ai + g log ^ + 0{a\ (62) 

where 

C' = hm I r ^ . : ^ + 



° -L-o I (27r)2 sin^ p^ + al + (r ^^{1 - cosp^)^ An 
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and 



2rE(l - cosp^) 



n 2 



In f l62|) . we discarded the leading constant that remains in the 0"^ ^ hmit, since 
it simply represents the uniform shift of vacuum energy. Note that the coefficient of 
a\ log 0"! agees with the correct one and the doubling phenomenon is not occured. Under 
(1591) . the linear term vanishes and thus the form of the effective potential agrees with 
that of the continuum limit, 



Vr 



2C2)al + ^\og^ + 0{a^). 



The dynamical mass of elementary fermion, is given by the solution of |^ 
From ([M]), we have the solution cr| = exp [47r(— + C'q 



(63) 
= 0. 



2g2 . -0 2C2)] = (amD)'-. 
Now, we turn to the issue that we have discussed through out this work. We try to 
capture the continuum scaling in expansions effective at large lattice spacings. In the 
case of the present model, we ask whether the large (Jl{= aa) expansion of the effective 
potential can be improved by the use of delta expansion or Borel transform. From ( l6Tl) . 



Vl at large cr^ read^j]] 



+ 



+ 



57r2 



r, , 4r 1 - / 28H\ 1 
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269rS 1 



1781r' 



6 



7r + 145r3 - 627r^ + 
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Fig. 11 shows the plot of V^af^ at r = 1/4, 1/2, 1 to the several orders in aj}. 
The dotted lines represent the asymptotic scaling of Vlct^^ (see f l63|) ) in the continuum 
limit. At r = 1, which is almost the standard choice in the literatures, the deviation of 
the exact function to the asymptotic behavior is small. For every sampled values of r, 
the asymptotic scaling begins around log ctl ~ —3. 

The delta expansion technique in the conventional scheme is not adequate in the 
present model. The difficulty originates from the existence of the linear divergence. The 
one-loop counter term cancels out the linear divergence involved in the integral, but the 
cancelation becomes perfect only when the integral is expanded in ctj^. When the integral 



II Though we have obtained I/ctl expansion from (ICT]) which holds only in the large N limit, we stress 
that the expansion is generally available through the hopping expansion. Here, we have used (j6ip only 
for the sake of efficiency. 
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Figure 11. Plots of the effective potential Vl/o^"^ with Wilson ferniion at g = f and 
r — i/4, f/2, f . The solid black lines represent the potential expanded in I/ul at 2nd, 
7th, 12th, • • • and 37th orders. The solid gray line represents the exact potential and 
the dotted line the asymptotic behavior ([55)1 in the continuum limit. At r = 1/4, the 
hopping series is effective to logCL ^ 0, while at r = 1, effective to logCL ^ 1. 



is expanded in l/a^, we cannot isolate the divergent piece and the conventional delta 
expansion makes the cancellation incomplete. The square type, on the otherhand, is 
useful to cope with the cancellation because all contributions are expanded to the same 
order of 5 in the scheme. Then, it is found that BT limit produces best simulation as in 
the previous models. Hence in the following, we report the result of applying BT limit 
to obtain the continuum limit of the effective potential. 

Now consider the Borel transform of V^o'l^- Using B[af^] = df^ jW. and 

B [-1 log al\ = (i log ai - + I) , (64) 

we obtain 

Fig. 12 shows the plot oiB^LOf^] for various values of r at orders 2, 7, 12, ■ ■ ■ , 37. When 
r = 1/10, the function oscillates at small a^. This property makes the quantitative 
analysis, for example the evaluation of mz), difficult. When r = 1/4, the approach 
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Figure 12. Plots of Borel transformed large ul expansion of Vlc^ and its asymptotic 
scaling at g = 1, B[VLal''] - - + 2^2 + ^(log 4 - 27b) (7b: Euler constant). 
From the first to the last, r = 1, 3/4, f /2, f /4, f /fO. Dotted line represents the 
asymptotic scaling and the solid lines \/(Jl series at 2nd, 7th, 12th, ••• and 37th 
orders. 



to the asymptotic scaling B\yLa~[^] ~ ~ + 2C*2 + 4^(log^ ~ '^Ie) is seen at 
12th order and the oscillation is weak. When r = 1/2, asymptotic scaling is captured 
around loga^ ~ —1.5 at 12th order and no oscillation is observed there. In general, the 
1/(Tl series becomes effective to smaller for smaller r but the expense is oscillatory 
behavior. As found from fl63p the leading logarithmic correction is independent of r 
but one cannot set r = from the outset. In other words, the l/o"/, expansion and the 
limit r ^ is not commutable. Oscillatory behavior at small r and small may be 
the signal of that non-commutablity. From quantitative point of view, most convenient 
value of the Wilson parameter would be a medium one around r ~ 1/2. 
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5. Conclusion 

We first introduced an alternative summation scheme called square scheme where the 
expansion in 5 is carried out to 5^ for all powers of M~^. Like the equalizer, the square 



shown that the delta expansion in square scheme leads to Borel transform in the limit 
M,N oo with M/N kept fixed. In the large J\f limit of A/" component anharmonic 
oscillator, this limit was shown to be natural by explicit application of square schemes. 

In anharmonic oscillators and non-linear a model, conventional and square schemes 
can be useful for the simulation of the continuum limit via 1/M expansion. In fact the 
accuracy of approximation obtained in BT limit is slightly worse than that obtained in 
conventional scheme employed in [3ll3]. However, it is important that in both truncation 
schemes, the double series in and 6 have been found to converge to the unique 

limit, the correct scaling, even in the 6^1 limit. In addition, flexibility and stability of 
the delta expansion under various truncation schemes was explicitly demonstrated. For 
example, from the connection to Borel transform, we can make rough estimation in what 
cases delta expansion is effective or not. One might guess that the expansion in 6 in large 
M series would not reflect the dilation around M = 0, if there exists a singularity on 
positive real axis of M. This is supported by the criterion of the Borel non-summability. 
Fortunately, lattice models allowing the continuum limit undergo second order phase 
transition and the transition occurs in the M — > limit and does not at M > 0. Hence, 
large M expansion belongs to Borel summable type and the delta expansion in various 
truncation schemes would help us to capture the continuum scaling in other complex 
lattice models. Thus, the mass in momentum space is a suitable parameter to express 
physical quantities. 

2d Gross-Neveu model is an example where the conventional truncation does not 
work and the square truncation, especially its BT limit, is necessary. This is because the 
linear divergence due to the Wilson term must be cancelled by the perturbative counter 
term while the original quantum correction needs to be expanded in l/a^. Using Borel 
transform, we found that the value of Wilson parameter around 1/2 is suitable for the 
quantitative study of the scaling and the dynamical mass evaluation. 
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